Effect of a thin optical Kerr medium on a Laguerre-Gaussian beam 
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Using a generalized Gaussian beam decomposition method we determine the propagation of a 
Laguerre-Gaussian (LG) beam that lias passed through a thin nonlinear optical Kerr medium. The 
orbital angular momentum per photon of the beam is found to be conserved while the component 
beams change. As an illustration of applications, we propose and demonstrate a z-scan experiment 
using an LGq beam and a dye-doped polymer film. 



Laguerre-Gaussian (LG) modes were discovered soon af- 
ter the invention of the laser—. A LG beam with angular 
mode number I and transverse radial mode number p can 
be written as: 




X exp 



X exp (2p+ 1^1 + l)tan ^ y — 

X exp exp {—ikz) , (1) 

where is the beam waist radius, — kLUQ/2 is the 
Rayleigh length, uj{z) = uioi^ + z'^/zq)^ is the beam ra- 
dius at z, and is an associated Laguerre polynomial. 
For a given LG beam mode, Zr (or ujq) alone is suffi- 
cient to characterize the relative amplitude and phase of 
the electric field of the beam. When multiple beams are 
involved, it is often necessary to specify the waist loca- 
tions of each beam. For convenience we use the notation 
C-LGp{r, (j>, z — Zw', Zr) to describe a LG beam unambigu- 
ously, where z^ is the waist location on the z axis and G 
is a complex constant that gives the amplitude and the 
initial phase. 

The fundamental Gaussian beam LGq has been the 
most commonly studied LG mode in both theory and 
experiment. However recently higher order LG beams, 
especially those with higher angular mode number I, are 
attracting more attention. When I > I, the LG beam 
possesses well defined orbital angular momentum of Ih 
per photon^,, zero on-axis intensity due to the phase sin- 
gularity, and helical phase fronts. These properties make 
high-order LG beams useful in many applications such 
as optical manipulation of microscopic particles^., quan- 
tum information processing^'^ , orbital angular momen- 
tum sorting^ and imaging!, nonlinear wave mixing^, op- 
tical vortices^ and vortex solitonai^, scatteringii, and 
interferencei^. In the present work, we develop a theory 
for the propagation of an LG beam after it passes through 
a thin nonlinear optical Kerr medium. The results of our 
calculation can be used for z-scan measurements using 
LGq beam and optical limiting geometries. 

The index of refraction of an optical Kerr medium de- 
pends on the intensity I of the incident beam and the 
nonlinear refractive index 712, 

71 = no + 772/, (2) 

where hq is the linear refractive index. Gonsider a nonlin- 
ear Kerr-type sample that is placed at the position z = Zs 



in a LG beam E{r, <f), z) ^ Eq ■ LG^°^ (r, (j), z; z^). Assume 
the sample is "thin" such that the change of the diameter 
of the beam due to either linear diffraction or nonlinear 
refraction can be neglected within the sample. Such two 
requirements are satisfied if d « Zr and d « Zr/A$, 
respectively, where d is the sample thickness, and A$ is 
the nonlinear phase distortionl4. The electric field E' at 
the exit surface of the sample can be expressed as^^'^^ 



E' (r, 0, z,) = E (r, 0, z,) e-¥e-»^*('-.^.-=), 
where a is the linear absorption coefRcient, and 
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where 



A$o = -ceo 77-0 712 1 -Bo I 
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(4) 



(5) 



is a constant proportional to the maximum nonlinear 
phase change A$maa:(-2^s) in the sample at the focus 
(zs = 0) . When the incident beam is a LGg beam, 

A$o = l^^„^ax{Zs = 0). 

We generalize the so-called "Gaussian decomposition" 
methodic to analyze the propagation of a beam that tra- 
verses a nonlinear sample. The exponential in Eq. ([3]) 
can be expanded in a Taylor series as 



-iA<E.(r,0,^,) ^ ^ (-7A$(r,0,Z3)) 
m=0 
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The complex electric field of the incident beam after it 
passes through the sample can be written as a summation 
of the electric fields of a series of LG beams of different 
modes as 



oo p. 



E' (r, 0, 2;) = ^ ^ Cp,mLG^p°^ (r, 0, z - z^,„; z™) , 



m— p— 



(7) 



Rayleigh length, respectively, of the corresponding beam 
mode and Cp^m is the amplitude and phase of the compo- 
nent beam. These parameters are determined by letting 
z = Zg in Eq. ([7]) and making use of Eqs. ^ and ([S]). 
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and Cp 
and 



(2m+l)(Z^ + l) 
Dp.m ■ -Pp.m where Z is defined a,s Z — z^/zr, 
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and and Dp^m are determined through 
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where x is an arbitrary real variable. Our theoretical 
results are embodied in Eq. ([7]). 

It is worth noting that all the component LG beams 
have the same angular mode number Iq as that of the in- 
cident beam, which reflects the conservation of the pho- 
ton's orbital angular momentum. Therefore the effect of 
the Kerr material on the incident LG beam is to gener- 
ate LG beams of different radial modes. These results are 
important in applications that leverage mode sensitivity. 

The expression in Eq. ([7]) is especially useful when the 
nonlinear phase distortion is small enough so that only a 
few terms in the summation are needed to make a good 
approximation. To illustrate, assume that the incident 
beam is a LGq beam and the maximum nonlinear phase 
distortion in the sample at position Z is 

|A<I>nl 

e • (1 + Z"^) 

It is sufficient to keep the first two terms in the Eq. ^ 
and neglect the higher order terms, yielding 

E' (r, 0, z) « FoflLGl (r, 0, z; Zr) 
2 

+ ^-^^0,1-^^0 (r, 4',z- Znii; Zri) 

Fi^iLG\{r,(j),z - Zwi;zri) . (9) 
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This outgoing electric field includes the generated LGq 
and LG\ beam. 

It is interesting to calculate the on-axis normalized Z- 
scan transmittance^'^'^'^ 

|2 



T(Z,A$o) = 



\E' (r ^ O,0,z ^ oo)Y 



(10) 



\E' (r ^ 0, z oo) |A<i.o=or 

which characterizes the on axis light power transmitted 
though a small aperture placed in the far field. Applying 
Eq. we find 

T Z, A$o = 1 + ^ A$o 

(1 + Z2)(9 + Z2)3 

.A$2. (11) 



(9 + Z2)^ 



The last term can be dropped if |A$o| << 1. 

The same procedures can be applied to LG beams of 
other modes. As a special case, when the incident beam 
is LGq, we get the same result as in Ref. [13, i.e. 



T(Z,A$o) = l 
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(1 + Z2)(9 + Z2) 
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(12) 



when |A(I>o| << 1- 

Result (IT^ was first reported as the theory support- 
ing the Z-scan technique, a highly sensitive technique to 
measure the optical nonlinearities using a LG^ beami^. 
In a Z-scan, one measures the trace of the normalized 
transniittance T as a function of the nonlinear sam- 
ple position z, from which A$o is determined by fit- 
ting Eq. (|12p or more practically, through the relation 
ATp_^ = 0.406|A<I>o| (derived from Eq. where 
ATp_^ is the difference between the peak (maximum) 
and the valley (minimum) of the trace (e.g., the dotted 
curve in Figure [21). A$o then gives n2 through Eq. 

Similarly, Eq. (fTTj) shows the relationship between T 
and A$o, suggesting that it is possible to do a Z-scan ex- 
periment using the LG\ beam to measure the nonlinear- 
ity of a thin samples. We demonstrate this on a disperse- 
red-1 doped poly(methyl methacrylate) (DRl/PMMA) 
thin sampleiS,. The experiment also tests the validity of 
our theory. 
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FIG. 1: Schematic diagram of the Z-scan experiment using 
a LGq beam. M: computer generated phase mask, L: lens, S: 
sample. A: aperture, and D: detector. 

Figure [T] shows a schematic diagram of the setup. A 
LGg beam of 647 nm wavelength from a krypton laser 
is converted into a LGq beam by a computer generated 
binary phase mask^^. The LGq beam is then focused 
by a lens to the desired Rayleigh length (2.4 cm). The 
power of the incident beam is 65 /iW. A 1.4 mm-thick 2% 
w/w DRl/PMMA sample is placed near the beam waist. 
The nonlinearity of DRl /PMMA is due to photo-induced 
trans-cis-trans isomerization of DRl molecules followed 
by reorientation in the direction perpendicular to the po- 
larization of the incident laser beam. It can be treated as 
an optical Kerr effect when the intensity is low and the 
exposure time is short enough to avoid saturation of the 
refractive index change. We use a shutter (not shown) 
to control the total exposure time (3 seconds). A small 
aperture is placed on axis of the beam in the far field. 
The power passed through the aperture is recorded by a 
detector as a function of sample position z as well as the 
exposure time t. The normalized transmittance is ob- 
tained by dividing the power recorded by the detector at 
time < = 3 s by that at time i = s. A fresh sample spot 
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is used for each exposure to avoid history effects. Non- 
Unear absorption is determined to be neghgible with an 
open aperture Z-scan. The "thin" sample approximation 
is satisfied under the above experimental conditions. 
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FIG. 2: Experimental (circles) and theoretical (solid curve) 
results of the Z-scan of a 1.4mm-thick 2% DRl/PMMA sam- 
ple using a LGl beam. Also shown are a LGq Z-scan trace 
(dotted curve) and a dashed line indicating the level of T = 1. 

Figure [2] shows the results. The circles are the exper- 
imental data. The solid curve is the best fit using Eq. 
pT|) with A$o as the parameter. The dotted curve is 
the best fit using Eq. (fT2| . included as a LGq Z-scan 
trace for comparison. The dashed line is to indicate the 
line where T = 1. 

The shape of the Z-scan curve using the LGq beam dif- 
fers from the traditional one using the LGq beam in that 
the former has an extra peak and valley (indicated in the 
figure by the arrows). The extra valley brings down the 
tail of the major peak below T = 1 and the extra peak 
brings up the tail of the major valley beyond T = 1. The 
experimental data clearly shows these features and fits 



our theory well. On the other hand, the tails in the LGq 
Z-scan trace never cross the T = 1 line and do not fit the 
data. Thus our theory provides the only good theory, 
supporting its validity. The appearance of the additional 
structure of our theory is a consequence of the LG\ com- 
ponent beam in Eq. ^ . Although the amplitude of these 
features is much smaller than the major ones, it has sig- 
nificance in other applications such as optical limiting, 
which will be discussed in a future publication. 

Just as the traditional LGq Z-scan, the LGj Z-scan 
provides sensitive measurements of optical nonlinearities. 
The value of A<i>o is fitted to be —1.1 (or A^maxiZ — 
0) = —0.4) using the data in Figure [21 With this we 
calculate n2 « —1.3 x 10"'* cm^/W, which is consistent 
with the value measured using the traditional LGq Z- 
scan. Besides using the fitting method, we can also es- 
timate A^„iax{Z — 0) using the difference between the 
major peak and valley, ATp_„ = 0.789\A^maxiZ = 0)|. 
An important difference between the two Z-scan mea- 
surements is that in the LG^ Z-scan experiment the de- 
tector is placed at the beam center where the intensity 
is the weakest due to the screw phase dislocation while 
in the LGq Z-scan experiment the center intensity is the 
strongest. As a result the former shows a much bigger 
deviation from the normal value if any phase or intensity 
distortion that destroys the symmetry of the beam profile 
is present. This suggests that the LGq Z-scan experiment 
is more sensitive to changes of the n2 of the sample. 

In conclusion, we have calculated the propagation of 
a Laguerre-Gaussian beam after it passes through a thin 
nonlinear optical Kerr medium. Applications, such as the 
Z-scan experiment using a LGq beam described here, has 
advantages over traditional measurements using a LGq 
beam. These results can be applied to many other thin 
film applications such as optical limiting and nonlinear 
beam interactions. 
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